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Abstract. We report variational and diffusion Quantum Monte Carlo ground- 
i , state energies of the three-dimensional electron gas using a model periodic Coulomb 

interaction and backflow corrections for N=54, 102, 178, and 226 electrons. We remove 
finite-size effects by extrapolation and we find lower energies than previously reported. 
, Using the Hcllman-Fcynman operator sampling method introduced in Phys. Rev. 

Lett. 99, 126406 (2007), we compute accurately, within the fixed-node approximation, 
the separate kinetic and interaction contributions to the total ground-state energy. 
The difference between the interaction energies obtained from the original Slater- 
determinant nodes and the backflow-displaccd nodes is found to be considerably larger 
^) , than the difference between the corresponding kinetic energies. 
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1. INTRODUCTION 

The homogeneous electron gas (HEG), which represents the simplest possible prototype 
of a many-fermion system, has been over the years a topic of intense research, as it 
often provides a good approximation for the description of valence electrons in simple 
metals and represents the basic ingredient for local and semilocal density-functional 
approximations. 

One of the first exhaustive calculations of the ground-state energy of an interacting 
three-dimensional (3D) HEG was performed by Ceperley PQ, using stochastic numerical 
methods. In this work, Ceperley used variational Monte Carlo (VMC) to obtain an 
upper bound to the ground-state energy. More accurate ground-state energies can be 
computed by using the more sophisticated diffusion Monte Carlo (DMC) approach, 
which projects out the true ground state of the system from a trial wave function [2]. 
However, this method yields a bosonic ground state even for a fermionic system which 
needs to have an antisymmetric ground state. In order to overcome this problem, the 
fixed-node (FN) approximation [3] has been applied; this approximation constrains the 
nodes of the ground-state wave function to those of a trial wave function. In spite of this 
constraint, the FN method has proven to be useful for the calculation of the ground- 
state energy and other electronic properties for atoms [HE], molecules [SIE], solids 019] 
and the two-dimensional (2D) electron gas pUl EU E2J- As an alternative to overcome 
the sign-problem for fermions, Ceperley and Alder [13] developed the so-called released- 
node DMC, which has the limitation that statistical fluctuations grow very rapidly and 
statistical noise can dominate the signal even before converging to the ground state. The 
size of the system that can be simulated is also limited in this method. Nonetheless, 
these released-node data have been widely used in the framework of density-functional 
calculations. 

With the improvement of computing capabilities and algorithms, new calculations 
have been attempted in order to improve the Ceperley-Alder DMC data. Ortiz and 
Ballone [H] extended these calculations to larger system sizes and polarized systems. 
Kwon et al. [15] introduced backflow and three-body correlations in the wave function, 
showing an improvement in the FN result beyond that given by Slater- Jastrow wave 
functions for a given finite system (see also [IE], [E] and [IE])- Both Ortiz and Ballone 
[13] and Kwon et al. [15] studied the size dependence and extrapolated their results 
to the thermodynamic limit; however, they assumed that the size dependence for VMC 
and DMC is the same, and they used VMC data to extrapolate the corresponding DMC 
calculations. 

Fixed-node DMC is known to yield the exact ground-state energy of a many-electron 
system for a given nodal structure. Nevertheless, the DMC expectation value of any 
operator that does not commute with the Hamiltonian differs from the exact value, the 
error being linear in the difference between the trial and the projected wave function. 
Recently, an effective method based on the Hellman-Feynman (HF) theorem was devised 
to calculate the exact expectation value of such an operator [19] as, for example, the 
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interaction energy. 

In this paper, we use the HF sampling introduced in [19] to report benchmark 
DMC calculations of the two (kinetic and interaction) separate contributions to the 
ground-state energy of a paramagnetic 3D electron gas with r s = 2^ We include 
backflow correlation effects and we demonstrate that previously extrapolated results 
are artificially lowered by assuming that the VMC and DMC size dependences coincide. 
However, we find that DMC size dependences are similar with or without backflow 
correlation effects. The use of the HF sampling allows to compute the exact (within 
the fixed node approximation) interaction energy in the framework of the modified 
periodic Coulomb (MPC) scheme, and we show that this interaction energy can also 
be obtained from the integration of the spherically averaged wave-vector dependent 
diagonal structure factor [20]. Furthermore, we find that the HF-sampling scheme works 
efficiently even for very large systems. 

Hartree atomic units (a.u.) are used throughout, i.e., h = \e\ = m e = 47re =l. The 
atomic unit of energy is e 2 /ao = 27.2 eV, ao being the Bohr radius. All the calculations 
presented in this work have been performed by using the casino code [21] . 

2. QUANTUM MONTE CARLO 

In VMC the ground-state energy, E VMC , is estimated as the expectation 
value of the Hamiltonian with an approximate trial wave function, ^t'- 
Evmc — (H)vmc — (^t\H\^/t)/ (^t\^t)- The integrals are evaluated by importance- 
sampled Monte Carlo integration. The trial wave function contains parameters, whose 
values are obtained from an optimization procedure formulated within VMC. There are 
no restrictions on the form of the trial wave function, and VMC does not suffer from 
a fermion sign problem. However, the choice of the approximate trial wave function is 
very important, as it directly determines the accuracy of the calculation. We have used 
VMC methods mainly to optimize the parameters involved in the trial wave functions 
by variance and energy minimization; our most accurate calculations, however, have 
been performed within DMC. 

In DMC the ground-state component of a trial wave function is projected 
out by evolving an ensemble of electronic configurations using the imaginary- 
time Schrodinger equation. The fermionic symmetry is maintained by the fixed- 
node approximation [3], in which the nodal surface of the DMC wave function is 
constrained to equal that of the trial wave function. The fixed-node DMC energy, 
E FN = (H) FN = (^ N \H\^ N )/{^ N \^ N ), is higher than the exact ground-state 
energy (H) = {^o\H\^/ ) / (^ol^o), and they become equal only when the fixed nodal 
surface is that of the exact ground state, \l/o- Apart from the fixed-node error, DMC 
yields the true ground-state energy independently of the form chosen for the trial wave 
function. The fixed-node error can be reduced by optimizing the nodal surfaces of the 

| r s is the dimensionless parameter r s = a/ao where ao is the Bohr radius and a is the radius of a 
sphere that encloses one electron on average. 
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trial wave function. 

2.1. The Slater- J astrow backflow trial wave function 

The standard Slater- Jastrow (SJ) wave function can be written as 

^sj(R) = e J ( R) ^ s (R), (1) 

where R is a 3iV-dimensional vector denoting the position rj of each electron. The nodes 
of ^sj(R) are defined by the Slater part of the wave function, \I/s(R)> which takes the 
form \l/s — -Df Ah At being a Slater determinant of single-particle orbitals of spin a. 
In the case of a HEG, these orbitals are Hartree-Fock solutions for the finite periodically 
repeated electron gas, which are simply plane waves. The number of electrons N has 
been chosen such that the ground state is a closed shell configuration, so the wave 
function can be chosen to be real and there is no degeneracy. The Jastrow correlation 
factor, e J ( R ), contains an electron-electron and a plane-wave term, as described in [22J. 
We did not include a symmetric three-electron Jastrow termj§| The Jastrow factor, being 
a positive definite function, keeps electrons away from each other and greatly improves 
wave functions in general, but it does not modify the nodal surface of the wave function. 

One way of reducing the FN error is to alter the nodes of the wave function by 
introducing backflow correlations [16J, thus replacing the coordinates R in the Slater 
part of the wave function by the collective coordinates X. The Slater- Jastrow backflow 
(SJB) trial wave function reads 

* SJB (R) = e J W* s (X). (2) 

The new coordinates for each electron are given by 

xi = ri +ei(R), (3) 

£i being the backflow displacement of particle i, which depends on the position of every 
electron in the system. Details of the specific form of the backflow function used for the 
HEG can be found in |16|. 



2.2. Hellman-Feynman sampling 

Given an arbitrary operator O, the fixed-node DMC method yields by construction the 
normalized expectation value (O)fn-dmc — (^t\0\^o N ) / (^tI^™) , which is not the 
true FN ground-state expectation value (0) FN = (^% N \d\ty£ N )/ (^% N \^% N ), unless the 
operator O commutes with the Hamiltonian, the leading term of this error being linear in 
the difference between ty? and ty™ . In conjuntion with VMC, this error can be reduced 
by one order by using the so-called extrapolated estimator, 2(0)fn_dmc ~ {0)vmc- in 
practice, extrapolated estimators work well when the trial wave function is very close to 
the ground state, but they can be untrustworthy when the trial wave function is poor. A 

§ We found that for the density considered in this work, the effect of the three-body term was not 
statistically significant in conjunction with backflow, which gives the best energy. This is in agreement 
with [16]. 
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correct sampling for these operators can be achieved by using future walking [23J 121] or 
reptation Monte-Carlo [25]; however, these methods aim at sampling ty™^™ instead 
of the regular DMC distribution [2], $r$Q 5V , so they are not straightforward additions 
to the DMC algorithm. 

An alternative to achieve a correct sampling of operators (diagonal in real space) 
that do not commute with the Hamiltonian has been reported recently [19]. The 
advantage of this so-called HF sampling is that it samples the usual DMC distribution, 
*&t^o N , and it is, therefore, straightforward to implement on a DMC algorithm. We 
give here a sketch of the method; a detailed derivation can be found in [19]. Given a 
Hamiltonian of the form H(a) — H + aO and the associated fixed- node ground-state 
energy E FN (a) = (H(a))pN, first-order perturbation theory for iff™ yields a fixed-node 
equivalent of the HF theorerrjjj] 

(6> FJV = ^M| . (4) 
da I a=o 

Direct application of the HF derivative to the regular DMC algorithm at timestep i 
gives: 

dEAa] 



O 



= Of-t{E{>X i -Et-X i ), (5) 

a=0 



da 

t is an auxiliary parameter and Of is the standard DMC estimator at time step i: 



() ! Yl u v ( 6 ) 

i 

where Uij is the total weight of walker j, and Ofj = O^t/^t with the trial wave 
function, ^t, evaluated for walker j at time step i. X t is the DMC estimator at time 

i 

step i of a new variable per operator X{j = | 0% j. The fixed-node DMC estimate, 

k=l 

which we call (O)fn-dmci is obtained averaging Equation fl6]) over all i. 

The correction term, AOf = —t ^Ef'Xi — Ef ■ XA , involves information which can 

be directly obtained from the DMC algorithm, such as the local energy, Efj = H^t/^t, 
and the new variable, Xi, which involves no more than an extra summation step during 
the sampling. The true FN estimate, which we call (O)hf, is obtained averaging Of 
over all i. 

Exponentially limiting the depth of the history in X^j allows one to considerably 
improve sampling and reduce statistical noise without reintroducing a significant bias 
[27]. 

3. Details of the calculations 

We have studied an unpolarized 3D HEG consisting of 54, 102, 178, and 226 electrons 
in a face-centered-cubic simulation cell subject to periodic boundary conditions. 

|| The Hellman-Fcynman theorem has also been applied previously in QMC for deriving accurate 
expectation values of observables [26] . 
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The parameters in the SJB trial wave functions were obtained by first minimizing 
the variance of the local energy [28| [29] and then minimizing the energy [30]. The specific 
forms of the Jastrow and backflow functions are described in [22] and [16]. The 2-body 
Jastrow term (U) and backflow (77) terms consist on power expansions in the electron- 
electron distance with expansion orders Nu=N v =8 and the parameters were allowed to 
depend on the spin parameters of the electron pairs. The cutoff lengths at which both 
U and 77 go smoothly to zero, Ljj and L v , were optimized, but they adjusted themselves 
to the maximum allowed values, i.e., the Wigner-Seitz radius. The plane- wave term in 
the Jastrow factor included 128 reciprocal-lattice vectors of the simulation cell. We used 
a sufficiently small time step (0.003 a.u.), to avoid finite-time-step errors, and a target 
population of 800 configurations in all our DMC calculations, making population-control 
bias negligible. 

Because our QMC calculations are performed using a finite simulation cell subject 
to periodic boundary conditions]^ the energy per particle is calculated at several system 
sizes and then the results are extrapolated to infinite system size. Finite-size effects in 
the kinetic energy are typically taken into account by noting that they are roughly 
proportional to the corresponding finite-size errors in the Hartree-Fock kinetic energy. 
Coulomb finite-size effects in the interaction energy of a HEG, which arise from the 
spurious interaction between an electron and the periodically repeated copies of its 
exchange- correlation (xc) hole, can be reduced either by adding the correction proposed 
by Chiesa et al. [33] to the usual Ewald energy or by using the MPC [3H [3U [351 EH] 
interaction. In this work, Coulomb finite-size effects were reduced by using the MPC 
scheme, where the Hartree energy is calculated with the Ewald interaction while the xc 
energy is calculated using 1/r within the minimum image convention, that is, reducing 
the interelectron distance into the Wigner-Seitz cell of the simulation cell. The MPC 
interaction was used for the branching factors and for computing energies in DMC. 
According to Drummond et al. [31] the Ewald interaction distorts less the xc hole 
when calculating branching factors in DMC. However MPC branching is faster and so 
potentially allows bigger systems. The trade-off is therefore between better convergence 
due to the xc hole being more accurate, or due to being able to go to bigger systems 
quicker. We have chosen the latter. 

4. Results and discussion 

Jf.,1. Ground-state total energy 

VMC and DMC ground-state total energies of a HEG of r s = 2 for iV=54, 102, 178, 
and 226 electrons, as obtained by using either SJ or SJB trial wave functions within the 
MPC scheme, are given in table [TJ These energies are also displayed by open symbols 

We used the Gamma point only. The recent analysis of Drummond et al. [31] showed that the 
T-point extrapolated SJ-DMC energies of a HEG agree within error bars with the results obtained by 
using the more sophisticated twist-averaged boundary conditions |32j . 
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in figure [T] with their corresponding error bars, which are less than 20% of the symbol 
size. Both VMC and DMC ground-state energies show the usual finite-size effects. For 
comparison, table [U also displays the DMC-SJ energies obtained by using the Ewald 
interaction. 



Table 1. VMC and DMC ground-state total energies, Evmc and Edmc-, as obtained 
with the use of SJ and SJB trial wave functions for a HEG of r s =2 for N— 54, 102, 
178, and 226 electrons and the MPC interaction. The column denoted with Ewald 
shows the DMC energies for the SJ wave function and the Ewald interaction. The 
entry N — oo corresponds to the extrapolated values obtained from Equation ((7|). 







SJ wave function 




SJB wave 


function 




MPC 


Ewald 


MPC 


N 


Evmc (Ha/e) 


Edmc (Ha/e) 


Edmc (Ha/e) 


Evmc (Ha/e) 


Edmc (Ha/e) 


54 


0.00847(2) 


0.00693(4) 


0.00426(2) 


0.00656(2) 


0.00579(2) 


102 


0.00328(1) 


0.00202(2) 


0.00081(2) 


0.001398(8) 


0.00077(2) 


178 


0.01057(1) 


0.00950(2) 


0.00892(2) 


0.008510(8) 


0.00801(1) 


226 


0.003490(8) 


0.00246(2) 


0.00202(2) 


0.001710(5) 


0.00123(1) 


oo 


0.003886(5) 


0.00301(1) 


0.00331(3) 


0.002021(4) 


0.001621(7) 



0.010 
0.008 



| 0.006 



0.004 
0.002 
0.000 




102 178 
Number of electrons (N) 



Figure 1. (Color online) Big open symbols: VMC and DMC ground-state total 
energies, Evmc and Edmc, from table [T] Horizontal dotted lines: The extrapolated 
values Eoo (also quoted in tabic [TJ, as obtained from Equation (J7J). Small filled symbols: 
The extrapolated values, En — b\ AThf(N) — 77, that we have obtained from Equation 
(JT]) for each V and fixed values of 61 and b^. 
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The issue of finite-size corrections to the kinetic and interaction contributions to 
the ground-state energy has been adressed by Ceperley and co-workers [TJ [37] Ell US] . 
They proposed separate extrapolation terms for the kinetic and interaction contributions 
using the following extrapolation formula: 



where E^, b±, and 6 2 are parameters to be fitted. ATn F (N) is the difference 
between the Hartree-Fock kinetic energies of the finite and infinite systems, 
ATn F (N) = Tft F (N) — Thf(oo), and the term b 2 /N accounts for the finite-size 
effects arising in the interaction energy. Recently it has been shown that the b 2 /N term 
also corrects for the neglect of long-range correlation effects in the kinetic energy |33j . 
Ortiz and Ballone [2] considered an extrapolation of the form 



with only 2 parameters to be fitted, and b . Av-hf(N) i n Equation (jSJ) is the 
difference between the Hartree-Fock exchange energies of the finite and infinite system, 
Avhf(N) = vhf(N) — Uhf(oo). El Equation ([TJ) and (jHJ) were found to yield similar 
results. We have found, however, that Equation (j7J yields in all cases better fits of 
the QMC data, so that all our extrapolations have been carried out by using Equation 
(JTjEI The horizonal dotted lines of figure [TJ display the extrapolated values E^ (also 
quoted in table [TJ that we have found from Equation flTJ). These extrapolated values 
indicate that the optimization of the nodes of the trial wave function (which is achieved 
by replacing the SJ trial wave function by the SJ-backflow trial wave function) lowers 
the ground-state energy considerably^ 2 mHa/e at the VMC level and 1.4 mHa/e at 
the DMC level; moreover, the optimized VMC-SJB ground-state energy happens to be 
1 mHa/e lower than its fixed-node counterpart DMC-SJ. 

Our extrapolated VMC and DMC ground-state total energies, E^, are compared in 
tableland figure [2] to the corresponding energies reported in [TJ [13], [HI [15] . These are: 
(i) the original VMC calculation of Ceperley [TJ , (ii) the released- node DMC calculation 
of Ceperley and Alder [13], (iii) the fixed-node VMC and DMC calculations of Ortiz and 
Ballone [H] and Kwon et al. [15], and (iv) the backflow DMC calculation of Kwon et 
al. [15] . We note that the DMC calculations reported by Ortiz and Ballone (OB-DMC- 
SJ)[H] and Kwon et al. (KCM-DMC-SJ and KCM-DMC-SJB) [15], were all carried 
out by first fitting within VMC the E^, bo, b\, and b 2 parameters entering Equation 
(J7J and (jHJ) and then using the VMC parameters b , bi, and b 2 to derive E^ for a given 
N from either Equation (j7j) or ([81: hence, in order to compare to the results reported 

+ Since the exchange hole entering the Hartrcc-Fock exchange energy is very long ranged compared to 
the exchange-correlation hole, the finite-size correction Avhf(N) entering Equation JHJ) might not be 
appropriate in an extrapolation scheme for QMC energies. 

* The adjusted R-squared values for the fits are: VMC-SJ: 0.999995, VMC-SJB: 0.999395, DMC-SJ: 
0.999004, and DMC-SJB: 0.999983. 

(j This corroborates the result reported in [16] to a small system of 54 electrons. 



E N = E oc + b 1 AT HF (V) + 



62 
V 



(7) 
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Table 2. Top: The extrapolated VMC and DMC ground-state total energies Eoo of 
table HI as obtained with the use of S J and SJB trial wave functions for a HEG of 
r s = 2. Middle: The VMC calculation of Cerperley (C-VMC) [I], the released-nodc 
DMC calculation of Ccpcrley and Alder (CA-DMC-RN) [13], the fixed-node VMC and 
DMC calculations of Ortiz and Ballonc (OB-VMC-SJ and OB-DM-SJ) [H] and Kwon 
et al. (KCM-VMC-SJ and KCM-DMC-SJ) [13], and the backflow DMC calculation of 
Kwon et al. (KCM-DMC-SJB) [T3J. Bottom: DMC ground-state total energies, but 
now obtained by using (as in [2] and [T3J) the VMC parameters b\ and 62 (which 
differ considerably from the corresponding DMC parameters) to derive Eqo for a given 
N from Equation (DMC-SJ* and DMC-SJB*). 





Eoo (Ha/e) 


Reference 


VMC-SJ 


0.003886(5) 


This work 


VMC-SJB 


0.002021(4) 


This work 


DMC-SJ 


0.00301(1) 


This work 


DMC-SJB 


0.001621(7) 


This work 


C-VMC 


0.002955 


ra 


CA-DMC-RN 


0.002055 


na 


OB-VMC-SJ 


0.0051(2) 


m 


OB-DMC-SJ 


0.0033(2) 


El 


KCM-VMC-SJ 


0.004096 


P3] a 


KCM-DMC-SJ 


0.002431 


P3] a 


KCM-DMC-SJB 


0.001812 


P3] a 


DMC-SJ* 


0.00266(1) 


This work 


DMC-SJB* 


0.001380(8) 


This work 



a These numbers for r s =2 were obtained by using an updated form of Equation (2) 
of [35], which was originally derived by fitting three values of r s (r s = 1, 5 and 10) 
of [15] and we have now fitted with the four available values of r s (r s = 1, 5, 10 and 
20) of p3] for each case (VMC-SJ, DMC-SJ and DMC-SJB). Differences between the 
energies derived from Equation (2) of [38] and the numbers reported here for r s =2 are 
within 7 x 10~ 5 Ha/e. 

by these authors we have performed additional DMC fits (DMC-SJ* and DMC-SJB* 
represented by dotted lines in figure [2]) by following this approximate extrapolation 
procedure. With this aim, we have calculated from the DMC-SJ (B) data and the 
b\ and 62 parameters of the VMC-SJ fit for each system size N, and then we have fitted 
a horizontal line to obtain the extrapolated DMC-SJ(B)* value reported in table El 

We find that (i) our DMC-SJ* calculation is very close to the corresponding 
calculation reported by Kwon et al. (KCM-DMC-SJ) [IS], and (ii) our DMC-SJB* 
calculation is considerably lower than the corresponding KCM-DMC-SJB calculation, 
which is a signature of the better quality of our SJB trial wave functions. We note that 
our calculations were performed using the MPC interaction, while in the calculations 
reported in the above references the Ewald interaction was used. Hence, we have 
repeated our DMC-SJ and DMC-SJB calculations using the usual 1/r Ewald interaction; 
these calculations (shown in tabled] for DMC-SJ) indicate that in both cases the Ewald 
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Figure 2. (Color online) Schematic representation of the extrapolated energies quoted 
in table H 



extrapolated energy is at most 0.3 mHa/e above the MPC extrapolated result, which 
reflects the fact that after extrapolation to N — > oo the usual Ewald energy yields 
fairly good results. We note, however, that the difference between the MPC and Ewald 
DMC-SJ energy reported in table [T]for each N is significantly smaller than the finite-size 
correction proposed in [33] (AV = u p /(4N), where u p is the plasmon energy). This is 
because the MPC and Ewald energies in tabled] were obtained using different interaction 
schemes in the branching part of the calculation, i. e., MPC and Ewald, respectively. We 
have performed additional calculations of MPC energies using the Ewald interaction in 
the branching factors, and we have found that only when the same interaction (Ewald) 
is used in the branching part of the calculation, the difference between the MPC and 
Ewald energies agrees with the correction proposed by Chiesa et. al [33] . 
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At this point, we focus our attention on a comparison between the DMC 
extrapolated values that we have obtained by (i) fitting within DMC the E^, bi, and 
62 parameters entering Equation ([7]) (DMC-SJ and DMC-SJB) and (ii) using the VMC 
parameters 61 and b 2 to derive for a given N from Equation ([7]) (DMC-SJ* and 
DMC-SJB*). Our calculations indicate that the size-dependences for VMC and DMC 
differ considerably^ indeed, the DMC-SJ* and DMC-SJB* extrapolated values are too 
low, i.e., the use of VMC data to extrapolate the corresponding DMC calculations 
yields artificially lowered extrapolations. Hence, the KCM-DMC-SJB ground-state 
energy reported by Kwon et al. [15] nearly coincides with our more accurate DMC 
extrapolation (DMC-SJB) as a result of two competing effects: the KCM-DMC-SJB 
extrapolated energy is (i) higher than our better optimized DMC-SJB* calculation and 
(ii) too low due to the assumption (in the extrapolation procedure) that the VMC and 
DMC size dependences coincide. 

Finally, we note that although the fitting parameters bi and 62 entering Equation 
(JTJ) are not transferable from VMC to DMC calculations they are indeed transferable 
from DMC-SJ to the more expensive DMC-SJB calculations: The error introduced by 
using the DMC-SJ parameters b\ and 62 to derive the DMC-SJB E^ for a given N from 
Equation ([7]) is found to be of no more than 0.04(1) mHa/e. 

4.2. Interaction energy 

The interaction energy, U, is a local operator (i.e., diagonal in real space) that does 
not commute with the Hamiltonian. Hence, for an accurate calculation of the true 
expectation value of U we have applied the HF-based method of [19] described in 
section 12.21 Figure [3] and table [3] show the results that we have obtained for the true 
(HF) fixed-node interaction energies Ufn = {U)hf of a HEG of r s = 2 for N = 54, 
102, 178, and 226 electrons, as obtained by using either SJ or SJB trial wave functions. 
Also shown are the fixed-node interaction energies Ufn-dmc — (U)fn-dmc that we have 
obtained by using the standard DMC estimator, which are subject to an error that is 
linear in the difference between ty? and ^q^. 

In order to obtain the extrapolated value Uoq from our finite-size calculations, we 
use the fitting equation 



for each set of data. The 1/N term does not aim at correcting the long-ranged errors, 
as we are using the MPC interaction. It turns out, however, that the residual effects 
including shell-filling and the distortion of the xc hole due to the finite size geometry 
also seem to be well described by a 1/N fit. The extrapolated values are displayed by the 
entry N = 00 of table [3] and the solid lines of figure [31 together with the result of using 



U N 




(9) 



ff Since the form and optimisation of the trial wave function entering VMC calculations is size- 
dependent, in general, there is no reason to expect VMC and DMC finite-size extrapolations to be 
the same. 



Benchmark Quantum Monte Carlo calculations 



12 



Table 3. Top: The fixed-node interaction energy of a HEG with r s =2 for N= 54, 
102, 178, and 226 electrons, as obtained by using the HF-based method of [H] (Ufn) 
and by using the standard DMC estimator (Ufn-dmc)] both SJ and SJB trial wave 
functions have been used. The entry N = oo corresponds to the extrapolated values 
obtained from Equation Bottom: The extrapolated estimator 2U fn-dmc — Uvmc 
for the infinite system. All energies are in Ha/e. 





SJ wave 


function 


SJB wave 


function 


N 


Ufn 


Ufn-dmc 


Ufn 


Ufn-dmc 


54 


-0.2993(1) 


-0.29839(9) 


-0.3009(1) 


-0.30050(9) 


102 


-0.29872(8) 


-0.2980(1) 


-0.3005(1) 


-0.2999(1) 


178 


-0.29782(7) 


-0.29741(6) 


-0.3001(1) 


-0.29958(9) 


226 


-0.2978(1) 


-0.29689(9) 


-0.2996(1) 


-0.2991(1) 


00 


-0.2973(1) 


-0.29679(8) 


-0.2994(1) 


-0.2990(1) 



2U D mc ~ U V mc (Ha/e) 

SJ (N = oo) -0.29764(8) 
S,m(N = oo) -0.2997(1) 



Table 4. The fixed-node extrapolated interaction energy ([/oo) for each system size 
as obtained from the fits of Equation (|9|) of the interaction energies of table [3] 





SJ wave 


function 


SJB wave 


function 


N 


Ufn 


Ufn-dmc 


Ufn 


Ufn-dmc 


54 


-0.2970(2) 


-0.2967(2) 


-0.2993(2) 


-0.2990(2) 


102 


-0.2975(1) 


-0.2971(1) 


-0.2997(1) 


-0.2991(1) 


178 


-0.29714(9) 


-0.29689(7) 


-0.2996(1) 


-0.2991(1) 


226 


-0.2973(1) 


-0.29648(9) 


-0.2992(1) 


-0.2988(1) 



the extrapolated estimator 2U fn-dmc — Uvmc for the infinite system (dotted and dashed- 
dotted lines). We see that while the true interaction energy Ufn is overestimated by 
the standard DMC estimator Ufn-dmc (the error being linear in the difference between 

an d *&q N ), it is underestimated by the extrapolated estimator 21) fn-dmc ~ Uvmc 
(the error this time being quadratic in the difference between ty? and ^q N )- On the 
other hand, we see that as in the case of the ground-state total energy the effect of 
backflow (included in the calculations labeled SJB) is to lower the interaction energy 
by a rigid shift of about 2 mHa/e. Finally, we note that error bars in figure [3] are 
approximately the size of the symbols (about 0.1 mHa/e), even for the largest systems 
under consideration. Table H] summarizes the extrapolated values for each system size, 
f^oo = Un — b/N, for all the cases shown in table [31 

In a recent paper [39], it was demonstrated that accurate calculations of the 
interaction contribution to the ground-state energy of an arbitrary many-electron system 
can be obtained from the knowledge of the spherically averaged wavevector-dependent 
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Figure 3. (Color online) Open symbols: The fixed-node interaction energies Ufn 
and Ufn—dmc quoted in table [3J Open symbols joined by solid horizontal lines: The 
extrapolated interaction energies obtained from Equation Black dotted and red 
dashcd-dotted lines: The extrapolated estimator 2Ufn-dmc — Uvmc for the infinite 
system, as obtained with the use of S J (black dotted line) and SJB (red dashed-dotted 
line) trial wave functions. 



diagonal structure factor Sf. as follows 

U= \J & -!]<**, (io) 
where S& is the spherical average of the diagonal structure factor in Fourier space: 

S k = 1 + drn(r) J duu 2 ^^-n xc (v,u). (11) 

Here, N is the particle number, n(r) is the electron density at r, and n xc (r,u) is the 
spherically averaged exchange-correlation hole density n xc (r, r 7 ) at r' around an electron 
at r. If one samples the structure factor using only correlations within the simulation 
cell, then Equation fllUI) represents the /c-resolved MPC interaction [20J. 

The structure factor can be computed using VMC [39], standard DMC, or the 
HF-based DMC [2D]. Since the spherically averaged structure factor is a diagonal 
function in real space, HF-based DMC calculations should yield the exact fixed-node 
S k . Figure 1 exhibits VMC, standard-DMC (FN-DMC) and HF-based DMC (FN) 
calculations of Sk for a HEG with r s = 2 and iV=102 electrons, as obtained with the 
use of SJ wave functions. Equation ([TP]) then yields the following VMC, standard- 
DMC, and HF-based DMC interaction energies: -0.29717(3) Ha/e, -0.2979(1) Ha/e, 
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and -0.29871(8) Ha/e, respectively, which agree with our calculated expectation 
values (U)vmc = -0.29719(3) Ha/e, (U) FN - DM c = -0.2980(1) Ha/e, and (U) HF = 
—0.29872(8) Ha/e (see also table E]), as expected. 




Figure 4. (Color online) VMC, standard-DMC, and HF-based DMC calculations of 
the spherically averaged structure factor Sk of a HEG with r s —2 and iV=102 electrons, 
as obtained with the use of SJ trial wave functions. 



4-3. Kinetic energy 

Table 5. Standard-DMC (top) and HF-based DMC (bottom) kinetic and interaction 
energies, as obtained in the thermodinamic limit (N — > oo) with the use of SJ and SJB 
wave functions. A = \SJ — SJB\ denotes the absolute value of the difference between 
the SJ and SJB calculations. 





t/oc (Ha/c) 


Too (Ha/e) 


FN-DMC SJ 


-0.29679(8) 


0.29980(8) 


FN-DMC SJB 


-0.2990(1) 


0.3006(1) 


AHSJ-SJB| 


22(1)10~ 4 


8(1)10- 4 


FN SJ 


-0.2973(1) 


0.3003(1) 


FN SJB 


-0.2994(1) 


0.3011(1) 


A=|SJ-SJB| 


22(1)10~ 4 


8(1)10- 4 



Assuming that the ground-state interaction and total energies have been correctly 
extrapolated, the kinetic energy can be obtained in the thermodynamic limit as 
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Tqo = Too ~ Uoq. Table shows standard-DMC and the more accurate HF-based DMC 
calculations of T^, as obtained from the corresponding standard-DMC and HF-based 
DMC calculations of Uoo (also quoted) with the use of either SJ or SJB trial wave 
functions. The absolute value of the difference between the SJ and SJB calculations 
(A = | S J — SJB | ) is also shown in this table. We note that this difference is considerably 
larger for the interaction energy, both in the case of the standard DMC approach and 
in the case of the more accurate HF-based DMC approach. This is an indication of 
the fixed-node error being smaller in the kinetic energy than in the interaction energy. 
Indeed, the correlation contribution to the kinetic energy is always smaller than the 
corresponding contribution to the interaction energy. 

5. Conclusions 

We have presented benchmark VMC and DMC ground-state energies of a 3D HEG with 
r s = 2 and iV=54, 102, 178, and 226 electrons, using an MPC interaction and backflow 
corrections. We have extrapolated our finite-size calculations to the thermodynamic 
limit, and we have found lower energies than previously reported, thus showing the 
good quality of our fixed-node trial wave functions. We have shown that previously 
extrapolated results are artificially lowered by assuming that the VMC and DMC size 
dependences (which we analyze independently) coincide. We have used the HF operator 
sampling method introduced in [19] to compute accurate values of the kinetic and 
interaction contributions to the ground-state energy. We also show that these values, as 
obtained with the use of the MPC interaction, coincide with the result one obtains from 
the spherically averaged structure factor. Our calculations indicate that our HF-based 
DMC approach yields very accurate results even for very large systems. Finally, we have 
found that the difference between the interaction energies that we obtain using either 
the original Slater-determinant nodes or the backflow-displaced nodes is considerably 
larger than the difference between the corresponding kinetic energies. A combination 
of (i) the fact that our Hellman-Feynman operator sampling method allows, within the 
fixed-node approximation, to calculate accurately the kinetic-energy contribution to the 
ground-state energy and (ii) the fact that the fixed-node error is smaller in the kinetic 
energy than in the interaction and total ground-state energy leads us to the conclusion 
that our kinetic energies should be of great use in the construction of accurate kinetic- 
energy functionals. 
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